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ABSTRACT 

This paper describes how parametric resonances associated with a galactic potential 
subjected to relatively low amplitude, strictly periodic time-dependent perturbations 
can be impacted by pseudo-random variations in the pulsation frequency, modeled as 
coloured noise. One aim thereby is to allow for the effects of a changing oscillation 
frequency as the density distribution associated with a galaxy evolves during violent 
relaxation. Another is to mimic the possible effects of internal substructures, satellite 
galaxies, and or a high density environment. The principal conclusion is that allowing 
for a variable frequency does not vitiate the effects of parametric resonance; and that, 
in at least some cases, such variations can increase the overall importance of parametric 
resonance associated with systematic pulsations. 

Subject headings: galaxies: evolution - galaxies: kinematics and dynamics - galaxies: 
structure 



1. Introduction and Motivation 

A real galaxy is never a completely isolated time- 
independent equilibrium. To a lesser or greater ex- 
tent it will always be subjected to a variety of per- 
turbations. These include (nearly) regular pertur- 
bations, reflecting, e.g., companion objects and/or 
(quasi-)normal modes triggered by close encounters 
with other galaxies {e.g., Vesperini & Weinberg 2000) 
as well as more irregular perturbations, reflecting 
e.g., internal substructures and/or larger numbers of 
nearby galaxies in a rich cluster. However, recent 
work {e.g., Kandrup, Vass & Sideris 2003) has demon- 
strated that periodic and/or near-periodic perturba- 
tions can trigger parametric resonances which make 
a significant fraction of the stars execute chaotic or- 
bits, leading to chaotic phase mixing and, oftentimes, 
readjustments in the bulk density distribution. 

Most work hitherto on the effects of such paramet- 
ric resonance has assumed a time-dependence charac- 
terised by a single fixed frequency. However, the fre- 
quency may not be exactly constant. Bulk changes in 
the density distribution during an epoch of violent re- 
laxation will cause changes in the bulk potential and, 
hence, the natural pulsation frequency. And similarly, 
for a galaxy near equilibrium, even if most of the time- 
dependence is associated with a single normal mode 
with a well-defined frequency, one may have to al- 
low for more irregular perturbations refiecting, e.g., 
nearby galaxies or internal substructures. 

In principle, these perturbations are determinis- 
tic. However, to the extent that they act in a 'near- 
random' fashion, one might wish to model them as 
coloured noise, i.e., random 'kicks' of finite duration. 
As emphasised in Siopis & Kandrup (2000), any noisy 
process can be viewed as a superposition of periodic 
influences combined with random phases. 

These considerations lead naturally to the issue 
considered in this paper, namely how parametric reso- 
nance associated with a single well-defined frequency 
can be impacted by perturbations inducing random 
fluctuations in the pulsation frequency. Do such fluc- 
tuations tend generically to enhance the effects of 
parametric resonance, or do they tend instead to 
weaken its eff'ects? In particular, are there circum- 
stances where even small random perturbations can 
dramatically increase the importance of parametric 
resonance and lead to potentially observable effects? 

One principal conclusion is that a variable fre- 
quency need not vitiate the effects of parametric res- 



onance; and, that, in some cases, such variations can 
make the effects of the resonance more pronounced. 
In particular, even for comparatively low amplitude 
perturbations, 'fuzzing out ' the frequency can lead to 
the displacement of appreciable quantities of stars 
and/ or gas to larger radii and, in some cases, to 
their ejection from the galaxy. This suggests that 
even comparatively weak interactions between galax- 
ies {e.g., the galactic 'harassment' of Moore, Lake & 
Katz [1998]) could contribute to the material often 
observed in intergalactic space in rich clusters, an ef- 
fect normally attributed {e.g., Bertin 2000) to galaxy- 
galaxy collisions and close encounters. This is an ef- 
fect likely to be missed in iV-body simulations, where 
the number of 'macroparticles' is small compared with 
the number of stars in a real galaxy. 

2. What Was Done 

As a simple example, one can consider orbits in a 
galaxy idealised as a spherically symmetric Dehnen 
(1993) model. Allowing for a time-dependent pertur- 
bation m{t), suppose that, in appropriate units, 

with — + + z"^ . Attention here will focus on 
7 = 0, 7 = 1, and 7 = 3/2, corresponding to galaxies 
with no cusp, a moderate cusp, and a steep cusp. 

Now suppose further that, in the absence of per- 
turbations, m{t) = 1, and, in a flrst approximation, 
allow for a periodic time-dependent perturbation 

m{t) = 1 + 5m{t) = 1 -|- mo sinwo*- (2) 

For suitable choices of mo and wq, such a perturbation 
will trigger a parametric resonance which can make 
orbits chaotic and, in some albeit not all such chaotic 
orbits, induce significant changes in energy. 

The objective then is to let the frequency to vary 
in a 'random' fashion which can be approximated as 
Gaussian coloured noise with a finite autocorrelation 
time. Assuming that Lxj{t) = wg + ^^{t)^ with 5u sam- 
pling a stationary Ornstein-Uhlenbeck process {e.g., 
Chandrasekhar 1943, Van Kampen 1981), the ran- 
dom process is characterised completely by the first 
two moments, 

{5uj{t)) = 

and 

{5uj{ti)5uj{t2)) = K{\tx-t2\), (3) 
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where ( . ) denotes a statistical average and the auto- 
correlation function 

K{\h - t2\) = A' exp(-|fi - talAc) (4) 

Here A represents the typical 'size' of the random 
component of uj and tc the time scale over which 
it changes appreciably. Alternatively, perhaps more 
physically, the noise process can be characterised by 
tc and a new quantity which corresponds to 

the mean value of 

The specific focus here was to imderstand how such 
random frequency variations change representative 
orbits evolved in an otherwise strictly time-periodic 
potential. The first task, therefore, was to explore 
how initial conditions corresponding, in the absence 
of perturbations, to purely regular, i.e., nonchaotic, 
orbits arc impacted by an exactly periodic driving. 
This involved determining the extent to which the 
values of radius r and energy E were affected by the 
perturbation. In the absence of the driving, energy 
E is of course conserved. However, in the presence 
of driving, this is no longer true and, for appropri- 
ate choices of mo and cjo, parametric resonance can 
induce large energy shifts. 

Having determined the effects of such periodic 
driving, the second, and principal, task was to se- 
lect 'interesting' values of unperturbed driving fre- 
quency Wo and explore how the dynamics is changed 
if Lo is perturbed by allowing for random noises with 
different values of mo, (li^wl), and tc- Most of the ex- 
periments assumed tc = 80. Noting, e.g., that, for 
the 7 = 1 Dchnen model, an intermediate energy or- 
bit with E = -0.4 has to ~ 10 (Merrritt & Fridman 
1996), this corresponds to a time longer than but still 
comparable to an orbital time scale. A more sys- 
tematic investigation of variable tc will be presented 
elsewhere. 

Since Su is a random variable, a single integration 
of an initial condition does not suffice . Rather, for 
each choice of initial condition, mo, ujq, and 
10000 different realisations were performed with the 
subsequent analysis focused on extracting statistical 
properties of the resulting orbits. The most impor- 
tant quantities to be extracted were the distributions 
oivmax and £^/i„, the largest value of r and final value 
of E assumed by each of the noisy orbits, which pro- 
vide information about the extent to which the noise- 
enhanced parametric resonance causes stars to be dis- 
placed to larger radii. (The quantity Vmax would seem 
a better quantity to track then r/m since, even if the 



(energy increases monotonically, the radial coordinate 
will oscillate; and it is quite possible that, at the end 
of the integration, r is near a 'relative' perigalacticon.) 
Attention focused separately on initial conditions cor- 
responding to (a) purely radial orbits and (b) orbits 
with a more isotropic distribution of velocities, ear- 
lier work (Terzic & Kandrup 2003) having shown that 
non-noisy parametric resonance can impact radial and 
nonradial orbits differently. 

This analysis is similar in spirit to recent work 
on noisy orbits in charged particle beams (Bohn & 
Sideris 2003), where the perturbations which one en- 
visions reflect, e.g., imperfections in the confining 
magnetic field and other details of an external en- 
vironment. However, the work described here is dif- 
ferent in at least two important respects. Firstly, the 
presence of a confining potential places an upper limit 
on the largest 'radius' that, in the absence of noise, 
a charged particle in the beam can achieve, which 
implies in turn that, for Su = 0, the response is a 
much smoother function of u>o than is the case for a 
galaxy. Secondly, and even more importantly, in the 
context of beam dynamics the principal issue of in- 
terest is in determining the largest values of r and 
E that any single beam particle can achieve, since 
the overriding concern is to prevent resonating parti- 
cles from striking the walls of the accelerator. In the 
context of galactic dynamics, the actual distribution 
of maximum radii and energies would seem more im- 
portant. Knowledge of these distributions can, e.g., 
facilitate predictions about expected changes in the 
density distribution and/or estimates of the number 
of stars that might be ejected into intergalactic space. 

3. What Was Found 

Consider first orbits subjected to strictly periodic 
driving. In this case, for fixed amplitude mo and 
autocorrelation time tc, the response to the time- 
dependence, as probed by Vmax and Efin, can ex- 
hibit a comparatively sensitive dependence on the 
frequency luq- Because the response reflects a res- 
onant coupling between (harmonics of) the driving 
frequency and the natural frequencies of the orbits, 
for sufficiently low and high luo the driving has a com- 
paratively minimal effect. In particular, the orbits are 
regular, i.e., nonchaotic, as was the case for the un- 
perturbed integrable model. For intermediate values, 
however, there are intervals where the orbits become 
chaotic and the response can be quite large. 
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For smaller amplitudes mo, the dependence of 
quantities like rmax on uq tends to be comparatively 
smooth, successive peaks in rmax{t^o), corresponding 
to strongly chaotic orbits, being separated by fixed in- 
tervals. However, for higher amplitudes the response 
can be more complex, especially for orbits in cuspy 
potentials. As noted elsewhere (Terzic & Kandrup 
2003), this is because Fourier spectra of cuspy orbits 
typically exhibit more structure than spectra of orbits 
in potentials with a smooth core. 

Examples of this behaviour for orbits in pulsed 
Dehnen potentials with 7 = and 7 = 3/2 are shown 
in Figure 1, which plots r,nax{^o) for initial condi- 
tions evolved for a time t = 512. For each poten- 
tial, results for two initial conditions are exhibited. 
Both correspond to an orbit with initial kinetic en- 
ergy T = —E/2. However, in one case the initial 
velocity had components Vx = Vy = Vz] in the other 
the velocity was selected to yield a radial orbit. In 
each case, the initial r = 0.25, this corresponding to 
a star situated near the center of the galaxy. For the 
cuspless potential, the responses for these two initial 
conditions are quite similar, but for the cuspy po- 
tential the radial orbit was often ejected to a much 
larger radius. This is consistent with Terzic & Kan- 
drup (2003), who found that, especially for cuspy po- 
tentials, radial orbits tended to be more susceptible 
to periodic driving. 

Now allow the frequency to vary by including col- 
ored noise in the computations, supposing that the 
unperturbed ivq is in one of the resonant regions. At 
least for 'reasonable' choices of parameter values, i.e., 
{\6(jj\) smaller than the width of the resonance, al- 
lowing the pulsation frequency to vary results on the 
average in orbits being displaced to larger radii. How- 
ever, when is very small, such variations have 
a comparatively minimal effect. If {\5w\)/u}o is very 
small, say < 10"'' or so, as a practical matter the fre- 
quency is 'nearly constant' in the sense that the varia- 
tions typically have no appreciable impact on orbital 
structure, so that nothing really changes. Alterna- 
tively, when becomes too large these variations 
can suppress the transport of stars to larger radii and 
less negative energies. If the variation is too large, the 
effective frequency with which the orbit is perturbed 
can be shifted out of the 'interesting' resonant range, 
at which point the driving no longer has a significant 
effect. The important point is that finite-sized vari- 
ations that wiggle the frequency within the resonant 
range tend generically to increase the mean rmax ■ 



Examples of this behaviour are exhibited in Figure 
2, which focuses on the same four initial conditions 
used to construct Figure 1, pulsed with frequency 
luq = 0.6 but allowing also for variations with vari- 
able (|(5a;|). In each case the data derive from 10000 
noisy integrations, each proceeding for a time t = 512. 
The solid curve represents the largest radial excursion 
experienced by any of the 10000 orbits. The dotted 
curve represents the mean value of rmax averaged over 
the different orbits. For very small and large values of 
(|5a;|), the frequency variations have only a compar- 
atively minimal effect but, for {\Su>\) ~ 10~^ — 10~^, 
the mean effect of this irregular time-dependence is 
substantially larger than for a purely periodic time- 
dependence. The range of values of {\Su}\) resulting in 
a significant effect is comparable for radial and nonra- 
dial orbits. However, especially for the cuspy poten- 
tial, within this range of values a varying frequency 
tends to have a larger effect on radial orbits. 

As is illustrated in Figure 3, a varying frequency 
can also increase the largest However, the mean 

value of Efin does not tend to grow significantly. 

Figure 2 demonstrates that the average effect of 
a variable frequency is to drive orbits to somewhat 
larger radii. However, it is also evident that some or- 
bits can be expelled to very large values of r: If the 
orbit and the noise are appropriately 'tuned', the reso- 
nant coupling can have an enormous effect. For exam- 
ple, a radial initial condition in the 7 = 3/2 Dehnen 
potential which, in the absence of a frequency varia- 
tions, is restricted to radii r < 1.2 can, in the presence 
of such variations, be ejected to r > 15! The obvious 
question, however, is how rare such events actually 
are. To address this issue, and to better understand 
the variations that a nonzero (|(5ti'|) can trigger, it is 
useful to examine f{rmax), the distribution of values 
of 

fmax for the 10000 different noisy realisations. 

Several examples of such distributions, generated 
for orbits in a 7 = Dehnen potential with = 0.6 
and mo = 0.05, are presented in Figure 4. Perhaps the 
most obvious point here is that these distributions can 
be quite complex. In general they are asymmetric, 
and they often exhibit multiple peaks, suggestive of a 
superposition of seemingly distinct populations. For 
the lowest amplitude variations exhibited in that Fig- 
ure, {\5u)\) « 0.0002, the mean values of rmax for the 
(?ns(^mbles arc very close to the r„iax computed in the 
absence of any frequency variations, but the distribu- 
tions f{rmax) still manifests significant structure. In 
particular, the distribution for the radial initial con- 
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dition exhibits two peaks, the existence of which is 
of definite statistical significance. For the intermedi- 
ate value {\duj\) « 0.06, the shape of the distribution 
f{r„iax) is relatively similar to the distributions in the 
upper panels. However, in this case most of the or- 
bits have managed to reach larger radii; only a few are 
restricted to values smaller than the r^ax associated 
with the driftless orbit. For the largest amplitude, 
{\5lo\) k, 0.25, even though the mean change in Vmax 
is relatively small, the distribution is quite broad; and 
it is clear that a majority of the orbits are actually 
restricted to smaller radii than was the driftless orbit. 

These effects may perhaps be better understood 
by focusing on changes in energy rather than radius. 
This is illustrated in Figure 5, which exhibits f{Efin), 
the distribution of final energies, computed for the 
same orbits used to generate Figure 4. Even visually, 
it is evident that, for the medium amplitude case, 
where the effect of the frequency variation is largest, 
f{Efin) is reasonably well fit by a Gaussian with a 
mean equal to the Efm associated with the zero vari- 
ation orbit. It thus appears that, at least as viewed in 
energy space, the variations result in strictly random 
changes in energy superimposed on the energy shift 
that arises in the absence of the variation. 

For the highest amplitude, the distributions are 
again Gaussian although the mean value is substan- 
tially smaller than the Efin for the driftless orbit. 
This decrease can likely be attributed to the fact that, 
in this case, is sufficiently large that, in many 

cases, the effective frequency acting on the orbit is 
outside the resonant region, thus yielding a weaker 
response overall. (As is evident from Figure 1 a and 
b, the width of the resonance near w = 0.6 is ^ 0.3.) 

For the lowest amplitude case, the situation is very 
different. For the isotropic initial condition, /(£'/j„) 
has an extended low energy tail; for the radial case, 
it is bimodal, indicative of two distinct populations. 
Given that the initial E = —0.24, it is evident that the 
total time-dependence (almost) never extracts energy 
from the orbit. However, the lower energy portion of 
the distribution corresponds to orbits which experi- 
ence only minimal changes. 

The key to understanding this behaviour is the fact 
that chaotic orbits in a strictly periodically pulsed 
Dehnen potential divide empirically into two 'types', 
namely 'sticky' chaotic orbits which exhibit little if 
any systematic changes in energy and 'wildly' chaotic 
orbits which do exhibit systematic changes (Terzic & 
Kandrup 2003). Analysis of Fourier spectra indicates 



that the 'stic;ky' orbits are 'locked' to a harmonic of 
the driving frequency, which restricts the extent to 
which they can drift in energy, whereas the 'wildly' 
chaotic are unconstrained. However, this distinction 
is not absolute. For example, an orbit that starts 
as 'sticky' will eventually become unstuck and start 
drifting through energy space in much the same sense 
that, in a time-independent Hamiltonian system, a 
chaotic orbit originally trapped near a regular island 
can diffuse through cantori to move into the chaotic 
sea. This analogy with time-independent Hamilto- 
nian systems suggests in turn that, just as noise can 
facilitate diffusion through cantori, corresponding to 
transitions between 'sticky' and 'wildly' chaotic be- 
haviour (Pogorelov & Kandrup 1999), the introduc- 
tion of weak noise into the driving frequency could 
facilitate transitions between orbits which do, and do 
not, exhibit systematic energy drifts. 

If the frequency shifts are large, it might seem un- 
likely for 'stickiness' to persist for long time intervals 
since this would require the orbital frequencies to vary 
in just the right way as to remain locked to the driv- 
ing. If, however, the frequency shifts are small, an 
orbit that started as sticky might be expected to re- 
main sticky for a relatively long time; and similarly, 
if the orbit starts as wildly chaotic, there is the pos- 
sibility that it will be 'kicked' into a frequency 'lock' 
and remain there for quite a while. 

That such behaviour is possible is illustrated in 
Figure 6, which exhibits E{t) and r{t) for an orbit 
evolved in a pulsed 7=1 Dehnen potential with ojq = 
0.6, mo = 0.01, and {\6u!\) = 0.003 for a total time 
t ~ 2048. Viewed over times t < 500 or so, the energy 
does not appear to exhibit any systematic changes 
and the radial motion appears comparatively regular, 
even though an examination of the Fourier spectra 
suggests that the orbit is not regular. It is, however, 
evident that, at a slightly later time, the energy begins 
to evolve in a more irregular fashion and that, for 
t > 1100 or so, the energy increases systematically. 

In any event, the behaviour exhibited in Figure 5 
for the medium amplitude case becomes progressively 
more common for larger amplitudes and cuspier po- 
tentials, where the resonances become broader and 
'sticky' chaotic orbits become less common. 

4. Conclusions and Interpretations 

Periodic driving tends generically to pump energy 
into the stars in a galaxy, thus displacing them to- 
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wards higher energies and larger radii, an effect re- 
sulting {e.g., Kandrup, Vass & Sideris 2003, Terzic & 
Kandrup 2003) from a resonant coupling between the 
driving frequency and the frequencies of the orbits, 
which also make the orbits chaotic. Even relatively 
weak driving with fractional amplitude ^0.05 can 
have significant effects within a few tens of dynam- 
ical times; the effects of larger amplitude can be large 
enough to account for violent relaxation (at least in 
principle). Given, however, that this shuffling in ener- 
gies involves a resonance, one might wonder to what 
degree this behaviour persists if the driving frequency 
is allowed to vary. 

One principal conclusion of this paper is that such 
variations do not vitiate the effects of such resonant 
couplings. Indeed, allowing for modest variations 
in frequency tends on average to increase the max- 
imum radii to which orbits are displaced. However, 
this does not imply that, on the average, such varia- 
tions result in more energy being pumped into the 
orbits. To the extent that the orbits have all be- 
come 'wildly' chaotic, i.e., that they have one or more 
positive Lyapunov exponents and that they are not 
'locked' to the (near-constant) driving frequency, al- 
lowing for different random frequency variations leads 
to a Gaussian distribution of maximum energies 
centered about the maximum energy attained by an 
orbit subjected to strictly periodic driving. The sys- 
tematic increase in the average r^ax arises because 
realistic (near-)equilibria have phase space distribu- 
tions which are monotonically decreasing functions of 
energy {df/dE < 0), so that a symmetric spread in 
energies occasions an increase in the number of larger 
radii stars at the expense of stars at smaller radii. 

More importantly, a random frequency can have 
a very large impact on at least a small number of 
stars. In particular, it is statistically probable that a 
few orbits will experience a noisy variable frequency 
which is 'just right' to give it a very large energy, 
thus displacing it to very large radii. (Alternatively, 
one might expect analogously that, for a fixed noisy 
variable frequency, there will be a small number of ini- 
tial conditions resulting in 'just right' chaotic orbits 
which will be impacted significantly.) In particular, 
the computations described above show that, in some 
cases, even a perturbation of fractional amplitude as 
weak as tuq 0.05 can cause a star to reach a ra- 
dius an order of magnitude larger than it would have 
reached if driven with a fixed frequency. 

This is an effect likely to be missed in numeri- 



cal simiilations of galaxy evolution. As illustrated 
in Figure 7 for Vmax, the distributions of values of 
Vmax and Efin decay exponentially for large values. 
This implies, e.g., that the largest value of Vmax com- 
puted for a collection of orbits should increase log- 
arithmically in N, an effect illustrated in Figure 8. 
Given, however, that few simulations of galaxies in- 
corporate more than ~ 10^ macroparticles, it would 
seem highly probable that such simulations could miss 
rare 'events' that arise at least occasionally in systems 
with > 10" stars. 

The computations described here suggest strongly 
that subjecting a galaxy to time-dependent perturba- 
tions which can be idealised as nearly periodic oscil- 
lations with variable frequency can have significant, 
potentially observable, effects. However, there remain 
a large number of lacuna, some of which are currently 
under investigation (Kandrup, Sideris & Bohn 2003). 
Most obvious, perhaps, is the question of how the re- 
sults derived here depend on the autocorrelation time 
tc- There are indications (Kandrup, Vass & Sideris 
2003) that the precise value of tc may not be all that 
important, but this remains to be checked. 

Also important is the issue of how coloured noise 
can facilitate transitions between 'sticky' and 'wildly' 
chaotic orbits. Earlier work (Terzic & Kandrup 2003) 
has shown that chaotic orbits subjected to purely 
periodic driving divide into two 'types', separated 
by 'entropy barriers', namely 'sticky' chaotic orbits 
which remain in the same phase space region without 
exhibiting systematic secular variations in energy and 
'wildly' chaotic orbits where changes in energy typi- 
cally grow diffusively. The obvious question, then, is 
to what extent the introduction of noise might ac- 
celerate transitions between orbit types. One knows, 
for example, that, in a time-independent Hamiltonian 
system, noise can facilitate diffusion through cantori 
or along the Arnold web (Lichtenberg & Wood 1989, 
Pogorelov & Kandrup 1999). 

Practical applications also require that the results 
described here be generalised to allow for unper- 
turbed potentials which are nonspherical, especially 
those admitting chaotic orbits. Earlier work (Kan- 
drup, Sideris, Terzic & Bohn 2003) has shown that 
time-independent potentials that already admit chaos 
can be significantly more susceptible to time-periodic 
perturbations, especially those characterised by rel- 
atively low frequencies (as might be appropriate for 
modeling an external environment); and it would not 
seem unreasonable to conjecture that frequency vari- 
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ations might also prove more important in such set- 
tings. 

Most importantly, however, it would seem imper- 
ative to consider the statistical effects of a noisy fre- 
quency on ensembles of initial conditions that con- 
stitute fair phase space samplings, so as to facilitate 
detailed predictions about the macroscopic properties 
of 'real' galaxies. In the first instance this could be 
done by considering e.g., microcanonical samplings of 
constant energy hypersurfaces in fixed potentials. Ul- 
timately, however, this problem can - and will - be 
addressed in the context of fully self-consistent simu- 
lations. 

Thanks to Court Bohn for useful discussions. HEK 
was supported in part by NSF AST-0307351. IVS was 
supported in part by Department of Education grant 
G1A62056. 
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Fig. 1. — (a) rmax, the maximum value assumed by 
the radial coordinate r for orbits generated from an 
isotropic initial condition with E — —0.24, evolved 
deterministically in a pulsed 7 = Dehnen potential 
with mo = 0.05 and variable frequency luq. (b) The 
same for a radial initial condition with E = —0.24. 
(c) The same isotropic initial condition evolved with 
mo = 0.01. (d) The same radial initial condition 
evolved with m = 0.01. In each case frequency was- 
sampled at intervals Sujo = 0.01. (e) - (h) The same as 
the preceding, now allowing for initial conditions with 
E — —0.55 evolved in a 7 = 3/2 Dehnen potential. 
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Fig. 2. — (a) The largest value (solid) and the mean 
value (dotted) of the maximum radial coordinate 
Tmax for noisy orbits generated from an isotropic ini- 
tial condition with E = —0.24, evolved in a 7 = 
Dehnen potential with luq = 0.6, nio — 0.05, and 
tc = 80. (b) The same for a radial initial condition 
with E — —0.24. (c) The same isotropic initial con- 
dition evolved with mo — 0.01. (d) The same radial 
initial condition evolved with m = 0.01. The dashed 
line corresponds to a deterministic orbit with Suj = 0. 
(e) - (h) The same as the preceding, now allowing 
for initial conditions with E = —0.55 evolved in a 
7 = 3/2 Dehnen potential. 
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Fig. 3. — The analogue of Figure 2 (a) - (d), now ex- 
hibiting the maximum energy Efin rather than rmax- 
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Fig. 4. — The distribution of maximal radial excur- 
sions, rmax for the two initial conditions of Fig. 1 
evolved in the 7 = Dehnen potential with luq = 0.6 
and mo = 0.05. The left panels are for the isotropic 
initial condition with (top to bottom) (|i5a;|) « 0.0002, 
0.006, and 0.25. The right panels are for the radial 
initial condition with the same values of The 
dotted vertical line exhibits rmax for the same initial 
condition subjected to purely periodic driving. 
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Fig. 5. The distribution of final energies Ejin for 
the same data analysed in the preceding Figure. The 
dotted vertical line represents Efin for the same initial 
condition subjected to purely periodic driving. 
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Fig. 6. — (a) The energy E{t) of a radial orbit with 
E{0) = —0.4 and y = z = evolved in a pulsed 7 = 1 
Dehnen potential with u>o = 0.6, mo = 0.01, and 
{\6lv\) = 0.003. (b) The coordinate x{t) for the same 
orbit. 
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Fig. 7. — (a) Exponential fall-off in the distribution 
fifmax) for an isotropic initial condition evolved in a 
pulsed 7 = Dehnen potential with u = 0.6, mo = 
0.05, and = 80, and {\6uj\) = 0.2. (b) The same for 
{\6lo\) = 0.007. 
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Fig. 8. — The largest value of Vmax for a radial initial 
condition evolved in a pulsed 7 = Dehnen potential 
with u = 2.8, m = 0.05, tc = 80, and {\Suj\} = 0.08, 
plotted as a function of the number N of noisy reali- 
sations. 
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